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Abstract 

We study a stochastic Schrodinger equation that generates a family of Gaussian dy- 
namical maps in one dimension permitting a detailed exam of two different definitions 
of non-Markovianity: one related to the explicit dependence of the generator on the 
starting time, the other to the non-divisibility of the time-evolution maps. The model 
shows instances where one has non-Markovianity in both senses and cases when one 
has Markovianity in the second sense but not in the first one. 



Recent theoretical and experimental advances have aroused a lot of interest in non- 
Markovian effects when quantum systems interact with an environment which cannot be 
considered at equilibrium pQ-pj)]. More specifically, consider a system S embedded in an en- 
vironment E, under the hypothesis of an initial factorized state, i.e., a density matrix of the 
form p® pe\ tracing away the environment degrees of freedom obtains an exact completely 
positive (CP) reduced dynamics for S that sends an initial state p at time t > into a state 
ptj at time t > t . This irreversible time-evolution is generated by an integro-differential 
equation of the form 

d t pt,t = duK tjU [p Uito ] , Ptoto =p, (1) 

J to 

where the operator kernel embodies the dependence on the past history of the system. The 
previous equation can be cast in the convolution-less form pU] 



dtPtM = U,t [pt,t ] , (2) 
where the presence of memory effects is now incorporated in the dependence of the generator 
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on the initial time t^. Because of this, the CP maps which solve (12"]) ,, 




Texp / dwJU )to , (3) 



1 t,t 

with T time-ordering, violate, in general, the (two-parameter) semigroup composition law, 
namely 

r t , tl °r tl ,^r tit0 , o<t <ti<*. (4) 

Indeed, if ~L u ,t = L u then (j3J) yields the equality in (J1J; vice versa, if in @ the equality 
holds, by taking the time derivative of both sides with respect to t one obtains L tjil = L ijto 
for all ti>t > 0. In [TO], the dependence of the generator L tjto on t and thus (|1J) is taken 
as a criterion of non-Markovianity. 

On the other hand, in p2] [T3] a different approach is considered whereby, given a one- 
parameter family of CP maps 7 t , t > 0, their non-Markovianity is related to non-divisibility, 
namely to the fact that no CP map A t >u , t > u > 0, exists that connects the maps 7$. In 
other words, the criterion of non-Markovianity becomes 

7i = A t , a o 7„ ^> A t>t0 not CP . (5) 

If a CP A i n existed, it would follow that certain CP monotone like the trace distance, the 
fidelity or the relative entropy should be decreasing: then, non-Markovianity is identified 
by the increase in time of such quantities which can also be taken as a measure of non- 
Markovianity. 

In order to study the two criteria of non-Markovianity, we consider a stochastic Schrodinger 
equation originally proposed as a non-Markovian mechanism for the wave function col- 
lapse [16]. Specifically, we take a particle in one dimension subjected to a time-dependent 
random Hamiltonian of the form (for sake of simplicity, in the following, vector and matrix 
multiplication will be understood) 

H? = H-w T (t)r, (6) 

where the Hamiltonian H is at most quadratic in position and momentum operators r T = 
(ri,r 2 ) = (q,p), while w T (t) = (wi(t),W2(t)) is a Gaussian noise vector with zero mean and 
2x2 correlation matrix D(t, s): 

=((w i (t)w j (s))), (7) 

j ij 

where ((■)) denotes the average over the noise. This latter matrix is real symmetric, Dij(t, s) = 
Dji(s,t), and of positive-definite type, that is 

D a(t*M > , V £(t«) e R 2 , (8) 
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for any choice of times {t a }™ =1 . For each realization of the noise, the Schrodinger equation 
(h=l) 



■ dm 

dt 



H-w T {t)r |C) 



(9) 



generates unitary maps U™ to on the system Hilbert space that send an initial vector state 
\ip) at time t = t into \ip™ to ) at time t. Averaging the projector \iptt ) i^ttol over the noise 
yields a density matrix 

ptM = im )m \)) ■ (io) 

In order to find U™ , one first goes to the interaction representation and sets: 



15.0) = ui to \^ to ), 

d|Co) 



dt 



^(t)f(t-to)l^l) , 



(11) 



where U t = exp(-iHt) and: 

r(t) = U}rU t = S t r , (12) 

St being a suitable symplectic matrix. For a given realization of the noise w(t), the solution 
is of the form l^uo) = ^t"t \^) where, a part for a pure phase, 



U™ = exp \ —i du w T (u) r(u — t c 



t 



w 

Mo 



(13) 
(14) 



By averaging over the noise, the corresponding density matrix ( TTOT) satisfies: 
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id t pt, 



to 



H, Pt 



to 



J2[^(Mt)Ko)Kto\)) 
3=1 



This stochastic Liouville equation can be turned into a standard master equation by means 
of the Furutsu-Novikov-Donsker relation [17J: 



w(s)X[w]}}= / du{(w{s)w{u))) 



,6R[w], 



(15) 



where -X"[io] is a functional of the noise, 5/5w(u) denotes the functional derivative with 
respect to the noise and R[w] is the density operator of the system. With R[w] = \?p™ to )(ipl 
one gets: 



w 

to I \ V y t,to I ' 



9tPt,t = L t ,t [p t . 



to J 



H, Pt 



to 



to J 



(16) 
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with: 
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^t,M = J2 <3y(Mo) {hpf, - \{r^p}) (17) 
C(t, t ) = / d« u) + Sl_ t D T (t, u)] . (18) 

If D(t,u) = 5(t — u) D (i.e., white noise) then one reduces to the Markovian Lindblad 
type dynamics with a time-independent positive Kossakowski matrix, namely C(t, to) = 
D [TH [19] . In the time-dependent case, in order that the maps T t j generated by L t>to be 
CP, the Kossakowski matrix C(t, to) need not to be positive, as we explicitly show in the 
following. We shall seek a solution of ffTB]) in the form 

Pt,t = T t , t0 [p] = J ^ G tjto (r) i2(r) W(5 t _ to r) , (19) 

where we have introduced the Weyl operators: 

W(r) = e irTn * = e i{qp -™ ) , (20) 

with r T = (q,p) G R 2 and Q — ^ ^ ^J, and i2(r) = Tr[p W(— r)] is related to the initial 
condition by: 

Pt„t =P = J ^R(r)W(r). 

Because the Hamiltonian H is at most quadratic and the matrix St in (fT2l) is symplectic, 
one finds: 

U t W(r)U} = W(S t r) . 
Direct insertion of (|T9|) into (|T6|) yields 



d t G t ,to(r) 



r T Sl t C(t,t )S t ^r CWr), 



whence G t ,t {r) = exp [-| r T gf(t, t ) r] with 
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flf(t,t ) = 2 / d M 5^ C( M ,t )5 u „ i0 (21) 

J to 

= [ du [ dvSl_ t0 D(u,v)S v - t0 . (22) 



Furthermore, since D(u,v) is of positive type, the matrix g(t,t ) is positive definite and 
Gt,t {f) a real Gaussian function; the solution I\ ito [p] can then be cast in a continuous 
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Kraus-Stinespring decomposition which guarantees the complete positivity of the maps T t to . 
Let G t>to (r) = J R2 d 2 x5(x - r) G t , t{) (x) with 

5(x-r) = -±- [ d 2 ye*^ n ^ . 
(2vr) 2 J R 2 

By inserting it into ( TT9l) and using W (x)W (r)W^ (x) = e~ ix nr W(r), one rewrites 

r tA) [p]=/ ^^(^^^(^p^t^)^ (23 ) 

Jr2 2tt 

with the Fourier transform 

Ft,M= [ ^jy Tnx G t , t0 (x), (24) 

also a real Gaussian, hence a positive function. 

Using (TIT)]) one can study the composition properties of the maps T tjto ; since: 

T^oIWpH J ^G t2 , tl (S tl ^r)G tuto (r)R(r)W(S t2 _ to r), 

in order to to satisfy the semigroup composition law r t2jtl o r tl to = T tl)to one should have 

G t2M (S tl „ to r) G tl>to (r) = G t2>to (r) . 
Using ( 1221) . one instead finds that 

/ / + / / dudv)[S T u _ to D(u,v)S v . t 
Jh Jt Y Jto Jt J \ 

I 2 du f 2 dv Sl_ t0 D(u, v)S v „ t0 . (25) 

Jto Jto 



This fact remains true even when D(s,u) = D{\s — u\) in which case from (122]) we have 

"t— to ft— to 



9{t,to) = / du dv S u D(u, v)S v 

Jo Jo 

and r t)to = T t _ to fi. 

Consider the master equation (fl6l) ; if to = its solutions p t fl = ^t,o[p] propagate the initial 
state p from t = to t > 0. Because of the above result, r ii0 ^ T tt t °^t ,o- However, setting 
t = in (fl6l) and searching a solution A tjto [p] in the form (fl9|) . one gets 

A M [P] = / ^L t , to (r)R(r)W(S t - to r) (26) 



where L tito (r) = exp { — ^r T £(t, t )r} with: 



£(t, tc 



/ duSl_ tQ C(u,0) S u - to 

J t 

pt pu 

du dv Sl_ tQ D(u,v) S v - 

J to JO 



t 



(27) 
(28) 



The function L tito (r) plays the role of Gt,t ( r ) hi (CGI to which it reduces when t = 0; that 
is A ii0 = r ii0 . Note however that, in contrast to g(t,t ) in (l2~Tj) . in £(t,to) one integrates 
C(w,0), not C(u,t ), from t to t. As a consequence, I\ i0 = A 4ji() o r 40)0 ; indeed, 

(l2r L t>t0 (S t0 r)L t0j0 (r)i2(r)W(5 t r) , 



K t ,t Q ° v to,M 
where now, unlike in (|25|) . 



2tt 



Sl£(t,t o )S to +£(t Q ,0) 



t PU 



to JO 



+ 



to pu 

Jo 



| d« dv D(u, v) S v 



du I dv S , D(u, v) S, 



£(t,0) . 

However, contrary to the maps A t)0 = T t0 which, as we have seen, are CP, the maps A t)to 
cannot be CP as this would imply [9] the positive definiteness of the matrix C(t,t ) in (TTTj) . 
In fact, the maps A t j are in general not even positive. 

All these various possibilities can be seen in a concrete example; consider a free particle 



of unit mass, H = p 2 /2, so that St 



given by 



1 t 
1 



, and a diagonal noise with correlation matrix 



D(t,u) 



- e -^- u \ fd q 
gL 



(29) 



2 \ w u,p 

First suppose the noise couples only to the position operator: d q — 1, d p — 0; then, from ([18 

1 _ e -7(*-*o) 



C(t,t ) 



e -7«-*o)[l +7 ( t _ to )]_l 
27 



e-7(*-*o)[i +7 (f_t )]-i ' 
27 







(30) 



has a negative eigenvalue for allt > t > 0. In spite of the non-positivity of the Kossakowski 
matrix in (|T8|) . the maps T t j in ( )23|) are nevertheless CP for all < to < 

We consider, as initial condition at to, a Gaussian state pa- with covariance matrix (CM) 
cr and zero first moments, Tr p a W{— r) = exp { — \r T (ft er f2 T ) r}. Using (TT9~j) . r tjto 

maps po- to the Gaussian state Tr r t ^ [pa-}W (r) = exp { — |r T f2 T cr t)to fir}, where cr t<0 = 
St-toVSj-to +fiit^o) with 

Sf(t, t ) = / dw fi T S£_ t C(u, t ) S u - t Q . (31) 

Jto 
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Instead, if the same initial condition is taken for the maps A t>to , the matrix g(t,to) is to be 
substituted by 

£(t, t ) = f du tl T S T u _ t C(u, 0) S u ^ t tl . (32) 

Jto 

If we choose <r = S to -t<^oSj _ t and expand a ti t = cro + &(t-,to) to first order about to, we 
have: 

<r tit0 ~<T + (t -t )n T c(*o,o)n, (33) 

where C(t ,0) is calculated from Eq. f[3"Ul) . Now, the second matrix at the l.h.s. is real 
symmetric and has one positive and one negative eigenvlaue, A > and — p < 0; let V be 
the symplectic, orthogonal matrix which diagonalizes it. Then, choosing an initial state with 
CM diagonal in the same basis, i.e., cr = Diag [u qq , cr pp ], such that cr + |fi > (positivity 
of the initial state), one gets: 

(Tt ' t °- V { a pp -^t-to)) V ' 

and a sufficiently small a pp would yield a non positive-definite CM a tt t , thus exhibiting the 
non-positivity of the map A t)to . The non-positive preserving character of A ttto is exposed by 
very specific states; on other states as, for instance, on all those of the form T to >0 [p] it acts 
perfectly well for A t:to o r to = I\ . In addition, starting from to = 0, A t = T t is CP. 

Therefore, in this case the master equation (Tl6i) generates a non-Markovian dynamics both 
according to the criterion since the generator L f) t depends on the initial time to and 
also according to the other criterion (jSJ). In fact, the family of maps IYo is non-divisible for 
A tt t Q is uniquely defined and non-positive. 

Since A ttto is not (completely) positive, certain quantities that exhibit monotonic behavior 
under CP maps fail to do so when evolving the system from time to to time t. One of such 
quantities is the fidelity [15] F{t) = J r (T t fi[p 1 ],T tt o[p2\) of two states p\ and p 2 evolving in 
time according to T^o- While J~{t) > J-"(0) for all t > 0, J 7 (to + t) may become smaller than 
J- (to) for some t, to > 0. This is showed in Fig. [TJ for two Gaussian states with zero first 
moments and "squeezed" CM. As one may expect, the effect disappears when 7 increases 
towards the Markovian limit. 

On the other hand, if in (Q, the noise affects the particle momentum only, namely if 
d q — 0, d p — 1, then, from (T2"5|) . 

C(t,t ) = (l-e-^) (° fj (34) 

is positive definite. It follows that the intertwining map At t t is CP, whence the family of 
maps T^o is divisible and Markovian according to the criterion (jSJ). However, it is non- 
Markovian according to the other criterion (j3J). Indeed, the generator resulting from (J9j) 
depends on the starting time t Q . 
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t 



Figure 1: Plot of the time evolution of the fidelity J- between two Gaussian states pk, k = 1,2, 
with zero first moments and CMs er^. = ^Diag [exp(2r^), exp(— 2r^)], with r\ = —r^ = 1.5, evolving 
under the map T^o for different values of 7. The inset refers to the time derivative of the fidelity. 
The non-monotonic behavior denotes non-Markovian evolution [15]; note that as 7 increases, T 
becomes monotonic. 



In conclusion, the analysis of above examples indicates that the criterion identifying non- 
Markovianity with the explicit dependence of the generator L tjto on the starting time t 
appears stronger than the criterion based on the non-divisibility of the maps T t Q . Indeed, 
on one hand, we have provided a case where the map T tt Q is divisible, yet the generator 
of Ft,t explicitly depends on the initial time to', on the other hand, a Markovian evolution 
according to the first criterion readily implies the semigroup composition law, i.e., (4) with 
the equality sign, hence divisibility of T t) o. Nevertheless, the non-divisibility criterion is the 
only one at disposal when one is presented just with the family of maps r tj o: in such a case, 
one may reconstruct the generator L tj o starting from to — 0, but, in general, no information 
is available on the full generator L t)to at t > 0. 
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